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Abstract 

Using a generalisation of the detailed balance for systems maintained out of equilibrium by contact with 2 reservoirs 
at unequal temperatures or at unequal densities, we recover the fluctuation theorem for the large deviation funtion 
of the current. For large diffusive systems, we show how the large deviation funtion of the current can be computed 
using a simple additivity principle. The validity of this additivity principle and the occurence of phase transitions 
are discussed in the framework of the macroscopic fluctuation theory. 

Resume 

Cumulants et grandes deviations du courant dans des etats stationnaires hors equilibre. 

En generalisant la relation de bilan detaille a des systemes maintenus hors equilibre par contact avec deux 
reservoirs a des temperatures ou a des densites differentes, nous retrouvons le theoreme de fluctuations pour la 
fonction de grandes deviations du courant. Pour de grands systemes diffusifs, nous montrons comment la fonction 
de grandes deviations du courant peut etre calculee simplement a l'aide d'un principe d'additivite. La validite 
de ce principe d'additivite et l'existence de transitions de phase sont discutees dans le cadre d'une theorie des 
fluctuations a Pechelle macroscopique. 
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1. Introduction 



A physical system in contact with two heat baths at unequal temperatures T a and Tt is one of the simplest 
situations for which one can observe a non-equilibrium steady state. At equilibrium, i.e. when the two heat baths 
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Figure 1. A system maintained in contact with two heat baths at unequal temperatures reaches in the long time limit a non-equi- 
librium steady state 

are at the same temperature (T a — Tb — T) , the probability P(C) of finding the system in a given microscopic 
configuration C is given by the usual Boltzmann-Gibbs weight 



P(C) = Z' 1 exp 



E(C) 



kT 



(1) 



where E{C) is the internal energy of the system in configuration C. Over the whole 20th century, studies in 
equilibrium statistical mechanics have been based on this expression or its microcanonical counterpart, and the 
great success of the theory was to show that (1) was the right starting point to explain the equilibrium properties 
of a large variety of physical systems (fluids, magnets, alloys, plasmas,....) and to understand all kinds of effects, 
in particular phase transitions and critical phenomena. A very simplifying aspect of (1) is that it depends neither 
on the precise nature of the coupling with the heat bath (at least when this coupling is weak) nor on the detailed 
dynamics of the system. 

As soon as the two temperatures T a and Tb are different [1], there is not such a simple expression [2,3] which 
generalizes (1) for the steady state weights P(C) of the microscopic configurations. In fact for a non-equilibrium 
system, the steady state measure P(C) depends in general on the precise description of the dynamics of the 
system, of the heat baths and on their couplings. So far the exact expression of these weights is known only for a 
few non- equilibrium models [4,5,6,7]. 

In addition to the steady state weights, one might be interested in the flow of energy through the system. For 
an interval of time t, one may consider the energy Q t , the energy transfered from the heat bath at temperature 
T a to the system. In the steady state, this energy fluctuates and one might try to predict its various cumulants 
(QDc or its large deviation function defined as 



Pro 



Qt 
t 



exp[-£F(j')] 



for large t 



(2) 



We refer to [8,9] for a full account on the large deviation theory. Note also that other definitions of the current 
distribution have been considered in [10,11]. 
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Figure 2. A typical shape of the large deviation function ^F(j) which vanishes at the typical current j. 
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The whole distribution of Qt and a fortiori its cumulants depend in principle on the initial configuration Cinitiai, 
on the final configuration Cfi na i and on the place where the flux of energy is measured. However if the internal 
energy of the system is bounded (maxc \ E{C)\ < oo), the cumulants of Qt (in the long time limit) and the large 
deviation function J-(j) do not depend on where the flow of energy is measured. In particular, if one measures the 
flux of energy between the system and the other heat bath, the large deviation function J-(j) is unchanged. Also 
if the system relaxes faster than the time t over which Qt is measured, the cumulants {Q") c divided by t and the 
large deviation function do not depend on the initial and final configurations Cinitiai, Cfinai- In fact in this 

case it is elementary to verify that J-(j) is convex, that is if < a < 1 

r(an + (1 - a)j 2 ) < aF{n) + (1 - a)T(j 2 ) (3) 

as the probability distribution Pro(Qt|Ci n itiai, C nna i) of Q t , given the initial and final configurations Cinitiai and 
C nna i, satisfies 

Pro(Q|Cinitial, Cfinai) = ^ ^ Pro(g| Cinitiai, C T )Pro(Q - q\C T , C't) > Pro( 9 | Cinitiai, C T )Pro(Q - q\C T , Cfinai) 
C T q 

which leads to (3) in the long time limit when r = at, q — jiat, Q — q = — a)t. The importance of convexity 
was understood in [12] (see Section 5 below). 

ft is sometimes easier to work with the generating function of Q t . For large t one has 

where (•} denotes the expectation over the dynamics and fx(\) is the Legendre transform of the large deviation 
function T 

/i(A) = max[jA - (5) 

3 

From the knowledge of /i(A), one can often determine the cumulants of Qt in the long time limit by 

(6) 



lim (QTh _ d n u(\) 



d\ n 

This relation is based on the assumption that the order of the limits t — > oo and A — > can be exchanged. One can 
show that these limits can be exchanged only for very few examples, although one believes that the assumption 
remains valid for general diffusive systems. There are however cases where these limits do not commute and for 
which the moments of the fluctuations cannot be deduced from the knowledge of the large deviation function [13]. 



2. Generalized detailed balance and the fluctuation theorem 



In principle determining the evolution of Q t requires the integration of the evolution equations of the system in 
presence of the heat baths. This is a difficult task, in particular because the heat baths are often described by an 
infinite number of degrees of freedom. Nevertheless, it can be shown in some cases that integrating the variables 
of the heat baths leads to effective reservoirs with stochastic noise. We refer to [14,15] and references therein for 
various ways of describing thermostats. 

Instead of considering mechanical systems, it is simpler to model the interactions with the heat baths by a 
stochastic term in the equations of motion of the system (like in a Langevin equation) . The microscopic dynamics 
becomes then stochastic. This means that the evolution is given by a Markov chain with transition matrix W(C' , C) 
which represents the rate at which the system jumps from a configuration C to a configuration C' (i.e. the 
probability that the system jumps from C to C' during an infinitesimal time interval dt is given by W(C' , C)dt). 

At equilibrium, one usually requires that the transition matrix satisfies detailed balance 

W{C',C) e — = W(C,C) e ~ (7) 

which ensures the time reversal symmetry of the microscopic dynamics. If one introduces q the energy transfered 
from the heat bath at temperature T to the system, and W q (C' ,C)dt, the probability that the system jumps 
during dt from C to C' by receiving an energy q from the heat bath, one can rewrite (7) 
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e* W q (C',C) = W- q (C,C") 



(8) 



If one accepts that detailed balance gives a good description of the coupling with a single heat bath at temper- 
ature T, one can wonder what would be the right way of describing the dynamics of a system coupled to several 
heat baths at unequal temperatures like in figure 1. When the system jumps from one configuration C to another 
configuration C' , energies q a ,qb,qc... are transfered from the heat baths at temperatures T a ,Tb,T c ... to the system. 
The straightforward generalization of (8) is 



q a , ib 

rp h Tpr ' 



.(C',C) = W- 



1a, -lb, -1c- 



XC,C") 



(9) 



For a system in contact with several reservoirs at temperature T a ,Tb,T c ..., this simply means, by comparing 
with (8), that the exchange of energy with the heat bath at temperature T a tend to equilibrate the system at 
temperature T a , the exchange with the heat bath at temperature T b tend to equilibrate the system at temperature 
T b and so on. 

The fluctuation theorem [16,17] can be easily recovered from the generalized detailed balance relation (9). To 
see this, one can compare the probability of a trajectory in phase space and its time reversal for a system in 
contact with two reservoirs. Similar approaches have been implemented for stochastic dynamics in [18,19,20,21]. 
A trajectory "Traj" is specified by a sequence of successive configurations Ci, ...Ck visited by the system, the 
times ti, ...tk spent in each of these configurations, and the energies q a ,i,qt,i transfered from the heat baths to the 
system when the system jumps from d to Cj+i. 



Pro(Traj) = dt k 



exp 



where r(C) = J2c E qa q W qa , qb (G" , C) and dt is the infinitesimal time interval over which jumps occur. 

For the trajectory " — Traj" obtained from "Traj" by time reversal, i.e. for which the system visits successively 
the configurations Ck, —Ci, exchanging the energies —q a ,i, —qt,i each time the system jumps from Cj+i to d, one 
has 



Pro(-Traj) = dt" 



HW- qati} - qbti (a,C i+1 ) 



exp 



E^(ft) 



One can see from the generalized detailed balance relation (9) that 



Pro(Traj) 
Pro {-Traj) 



= exp 



\ - q a ,i Qb,i 
^ T. T b 



— exp 



(a) 



T b 



(10) 



where Q^ a ' = q a ,i and Qf = q b ,i are the total energies transfered from the heat baths a and b to the system 



-,(6) 



during time t. If the internal energy of the system is bounded, energy conservation implies that \Qt^ + Q K t'\ < E, 
and one gets 



1W1 



exp 



Q 



(a) 



1 

n 



i 

Ta 



E_ 



Pro(Traj) 
Pro(— Traj 



< exp 



1 

n 



i 



(11) 



If P(C) is the steady state probability of configuration C, the probability that Q t = is the total energy 
transfered from the heat bath a to the system is given by 

Pro(Q t )= EE E ^(Cinitiai) Pro(TTOj(G nitia i,C' fina i,Qt)) 

^initial 

where the sums are over all initial configurations Cmitiai, final configurations Cfl na i and all trajectories Traj(Cinitiai, Cfi na i, Qt) 
starting in configuration Cinitiai, ending in configuration Cfinai with a total transfer of energy Qt. Now as 



Pro(-Q 4 )= EE E ^(Cflnai) Pro(Traj(C flna i,Ci ntia i,-Q t )) 

Cinitiai Cflnai Traj (Cf ina i ,C init i a i , - Qt ) 
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one can see that, if for any pair of configurations the ratio of their steady state weights remains bounded 
VC,C 0<.4< <B<oo , 



one has because of (11) that 
A 



s exp 



^ Pro(Qt) „ B 

< p7oF^) < 1 cxp 



E_ / 1 1_ 

T b +Qt \T b T a 



Taking the log and then the long time limit (2) leads to the fluctuation theorem 

m - H-j) = -3 (Jr - -r) (12) 

which states that the difference T{j) — is linear in j with a universal slope related to the difference of the 

inverse temperatures. 

We see that in the framework of stochastic dynamics, the fluctuation theorem is an elementary consequence of 
the generalized detailed balance relation (9) satisfied by the dynamics and of the assumptions that the energy is 
bounded (see [22,23,24] for examples where the energy is not bounded in which case the fluctuation theorem has 
to be modified) and the fact that the time t is much longer than the relaxation times in the system. In terms of 
the Legendre transform (4,5) the fluctuation theorem becomes 

M(A)=^(-A+^-^) (13) 
Remarks: 

(i) In the limit of small T a — T b (i.e. close to equilibrium), one can recover from (13) the fluctuation-dissipation 
relation between the variance of the current at equilibrium 

^jt -> 5- for T a = T b (14) 
and the response to a small temperature gradient 

{T a - T b )D for T a - T b small (15) 
In fact from these definitions of a and D, one has 

M(A) = {T a - T b )D\ + a -\ 2 + O (A 3 , \ 2 {T a - T b ),\{T a - T b f) (16) 

and for this expression to satisfy the fluctuation theorem (13), the coefficients a and D have to satisfy 

cr = ITlD (17) 

which is the usual Einstein fluctuation-dissipation relation between the response coefficient D and the 
fluctuation coefficient a. Note that in general both D and a depend on the temperature T a . 

(ii) One can easily extend the generalized detailed balance (9) and the fluctuation theorem (12,13) to other 
types of currents. For example, in the case of a current of particles, (9) becomes 

za q * z b - q » W qa , qb {C',C) = W- qa ,- n {C,C) (18) 

where z a and z b are the fugacities associated to the reservoirs of particles and q a and q b are the numbers of 
particles transfered from the reservoirs while the system jumps from configuration C to configuration C' . 
The fluctuation theorem (12,13) becomes then 

T{j) - = j[logZf, - logz a ] and /i(A) = n (-A + log z b - log z a ) (19) 

Close to equilibrium, if one defines as in (14,15), the fluctuation and the response coefficients for a system 
in contact with two reservoirs 

(O 2 ) (O ) 

> a for p a = p b and — j > {p a — p b )D for p a — p b small (20) 
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where D and a are now functions of the density p a . One can show, by expanding in powers of A and of 
z a — Zi as in (16) that 

a = 2D-^— = 2DTp 2 n (21) 
a log 2 

where k = p~ 1 dp/dp is the compressibility. (To see why the compressibility appears, one can write logZ = 
—W/T = —Vf(N/V)/T where F is the free energy and / the free energy per unit volume; one uses the facts 
that Tlogz = dW/dN — f'(p) and that p — —dW/dV = pf'(p) — f(p); then one can see that dp/dlogz — 

T/r(p) = T P dp/dp). 

(iii) Another easy extension is to consider systems with several types of currents (for example a current of 
particles and a current of energy, or several types of particles, or systems in contact with more than 
two reservoirs). The extension of the fluctuation theorem to these cases allows one to recover Onsager's 
reciprocity relations in the close-to-equilibrium limit [25,20]. 

(iv) The fluctuation theorem is usually formulated in terms of entropy production [26,27,21] as, in the steady 
state, the entropy of the system remains stationary whereas a current j of energy from the heat bath at 
temperature T a into a heat bath at temperature Tb gives a rate j(l/Tb — 1/T a ) of increase of entropy. 

An example: the symmetric simple exclusion process (SSEP) [28,29,30,31] 
There are only few examples of non-equilibrium steady states for which the cumulants (Q") c or the large deviation 
function J-(j) of the current can be calculated [32,33,34]. One of the simplest cases is the symmetric simple 
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Figure 3. The symmetric simple exclusion process 

exclusion process shown in figure 3. 

The model is defined as a one dimensional lattice of L sites with open boundaries, each site being either 
occupied by a single particle or empty. During every infinitesimal time interval dt, each particle has a probability 
dt of jumping to the left if the neighboring site on its left is empty, dt of jumping to the right if the neighboring 
site on its right is empty. At the two boundaries the dynamics is modified to mimic the coupling with reservoirs 
of particles: at the left boundary, during each time interval dt, a particle is injected on site 1 with probability adt 
(if this site is empty) and a particle is removed from site 1 with probability ^dt (if this site is occupied). Similarly 
on site L, particles are injected at rate 5 and removed at rate f3. (Note that one could consider that in the SSEP 
the particles represent quanta of energy and all the properties could be interpreted in terms of heat transport.) 

From the definition of the model, it is immediate to see that the dynamics satisfies the generalized detailed 
balance relation (18) with 

a 8 

Za = ; Zb = — 

7 P 

If n is a binary variable which indicates whether site i is occupied (t; = 1) or empty (t; = 0), it is easy to calculate 
the steady state profile [32] 

L + 1 + a + b 

where 

Pa = — ; — , pb = -5— j and a = — ■ — , b = — — . 

a+7 p + o a + 7 p+o 

Clearly, in the expression (22) of the profile, p a and pb represent the densities in the reservoirs at the two ends of 
the chain. The calculation of the first cumulants can be done either directly or by a perturbation theory in A by 
using the generating function [32]. 

lim iQlL = Pa Z Pb 
i^oo t L + a + b — 1 
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(Qt)c l,„. n n, „x , a(g-l)(2a-l) + b(ft-l)(2b-l)-Li(Li-l)(2Li-l) 2 
^ — = T^ pa + Pb ~ 2papb > + 3L?(L!-1) (Pa ~ Pb) 



lim 



where Li=L + a + 6— 1. For large L, the first four cumulants are given by 



lim 

t — *-oo 


(Qt) „ 
t 


. P a ~ 

L 




lim 

t — s-oo 


(Q 2 t)c 
t 


~ 1 
~ I 


pa 


lim 

t — *-oo 


(Q 3 th 
t 


~ 1 
~ I 


pa 


lim 

t — *-oo 


{Qih 
t 


~ 1 
~ I 


Pa 



2p 2 a + 2p aPb + 2pp 



Wpl + 12p 2 aPb - \2p a p\ - IQpl 



vro r "' 15 

lApj + 2p a p b + Up 2 b 32p 3 a + 8p 2 a p b + 8p a pi + 32p 3 b 
3 5 

96p 4 a + Mp 3 aPb - AQplpl + QApapl + 96pi 



3.1 



(23) 




Figure 4. The fourth cumulant versus for p a — 1. The thin lines represent the fourth cumulant obtained from exact calculations 
of /ii(A) for system sizes L — 5, 9, 13, 17, whereas the thick line represents expression (23) valid in the limit L — > oo. 



3. The additivity principle 



One can formulate a conjecture, the additivity principle [35], based on a simple physical picture, which allows one 
to determine all the cumulants and the large deviation function J-(j) for more general one dimensional diffusive 
systems. Applied to the SSEP, this leads to the same expression of the cumulants (23) and provide a way of 
calculating all the higher cumulants. Here we will limit the discussion to non-equilibrium steady states of systems 
in contact with two reservoirs of particles. As shown below everything can be easily generalized to systems in 
contact with two heat baths. For a system of length L + V in contact with two reservoirs of particles at densities 
p a and pi,, the probability of observing, during a long time t, an integrated current Qt = jt has the following form 
(2) 



Pro L+L / (j,p a ,Pb) ~ e~ 



(24) 



The idea of the additivity principle is to relate the large deviation function F L+L i(j, p a , p b ) of the current to the 
large deviation functions of subsystems of lengths L and L' by writing that for large t 



Pro 



L+L- 



(j, Pa, Pb) ~ max [Pro L (J, p a , r) x Pro L / (J, r, p b ) 



(25) 



This means that the probability of transporting a current j over a distance L + L' between two reservoirs at 
densities p a and p b is the same (up to boundary effects which give for large L subleading contributions) as the 
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probability of transporting the same current j over a distance L between two reservoirs at densities p a and r times 
the probability of transporting the current j over a distance V between two reservoirs at densities r and p b . One 
can then argue that one should choose for r the density which makes this probability maximum. From (25) one 
gets the following additivity property of the large deviation function 

Fl+l' (j, Pa,p b ) = min [F L (j, p a ,r) + F v (j, r, p b )] . (26) 

r 

By repeating this procedure, one gets that 



i=0 



F L {j,Pa,p b )= min {JZFi&r^n+i)} (27) 



where k — L/l, r — p a and r k — p b . 




1 X 

Figure 5. The dashed line represents the steady state profile. The density changes in the bulk to facilitate the deviation of the 
current. 

For large L and k (with L/k still very large), if one considers current fluctuations of order 1/L, it is advantageous 
to minimize (27) to make the differences n — n+i small. As the current j is also small one can consider that each 
piece of length I is close to equilibrium and has Gaussian fluctuations at the leading order 

r • _ -P(n)(n-n+l) i2 

F ( 0-,r i; r i+1 ) ~ y. — ^ (28) 

where the parameters D and a are defined for a system of length I in contact with two reservoirs at densities p a 
and pb by 

-> -y- for p a = p b (29) 

^-^(Pa-Pb) for p a - p b small (30) 

These are the same parameters as in (14,15,20) up to a factor I. (In the definitions (29,30), one should take first 

the t — > oo limit, i.e. cy{p a ) = limi^oo linit^oo ). 
If for large k, the density n varies slowly with i 



for some smooth density p(x) (see Figure 5) then combining (27) and (28), we get 

F L (j,p a ,p b ) = T n g =mm - J dx (31) 



with p(0) = p a and p(l) = pb. Note that (28) is a local equilibrium assumption, i.e. that both the current j and 
the difference n — r»+i are small. Therefore one cannot expect (31) to be valid when the current deviation j is 
not of order 1/L. 
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The average profile p(x) is the one which makes vanish the large deviation function Fl ■ It therefore satisfies 

Lj + D(p(x))p(xy = 
where the most likely current j is fixed by the boundary conditions 

Pa 

~3 = \\ D (p) d P- 



4. The large deviation function obtained from the additivity principle 

4.1. The optimal profile 

The profile po(x) which optimizes (31) satisfies 

d L 2 f _ , 2 d D 2 (po(x))) _ „ D 2 (po(x))) 
dp2a(p Q (x)) m 1 dp 2a(p (x)) m ' 2a(p (x)) 

If one multiplies this expression by p'o(x), one can integrate once. Finally one gets that the optimal profile satisfies 
(Lj) 2 (l + 2Ka(p (x))) 

M = Wm) (32) 

where the integration constant K is fixed by the boundary conditions po(0) = p a and po(l) = Pb- 

Suppose that p a > p b and that the deviations are not too large so that the optimal profile remains monotone, 
i.e. 



, Tj y/l + 2Ko(po(x)) 

p {x) = -Lj — (33) 

D(po(x)) 



one can rewrite (31) as 

we ^ ■ f \ 1 + K(J (P) 
F L (j,p a , Pb )= 3 J [ [1 + 2Ka{p)Y/ , -\ 



'^dp (34) 



Pb 

where the constant K is fixed from (33) by the boundary condition (p(0) = p a and p(l) — pt), i.e. 

"-/ pAp" * 1351 

Pb 

The optimal profile (33) remains unchanged when j — » — j (simply the sign of [1 + 2Ka(p)] L ^ 2 is changed) in 
(34,35) and one gets that 

Pa 

F L {j)-F L {-j) = -2j J ^ dp (36) 

Pb 

which is the fluctuation theorem (19). In fact already in (31) it was clear by expanding the square that the optimal 
po(x) does not depend on the sign of j and that (36) had to be satisfied. 

The physical meaning of the optimal profile po (32) is that adopting this profile is the easiest way to flow through 
the system an atypical current j. The large deviation functional (31) shows that the optimal density profile p 
and the current deviation j are coupled in a non trivial way. One can think of the system as a pipe with diameter 
<j(p) depending on the local density. The easiest way to increase the particle current is to adjust the size of the 
pipe a(p) and therefore the local density, in order to facilitate the flow of particles. In the example of the SSEP 
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with reservoirs at equal densities p — p a — p b , the variance of the current a(p) = 2p(l — p) is maximum at density 
p = 1/2. When p a = p b < 1/2, it is favorable to have in the bulk a density po(x) > p a in order to facilitate the flow 
of particles and the optimal way of doing it is by choosing the profile po(x) which satisfies (32). If p a = p b — 1/2, 
then the optimal density profile remains flat for any current deviation and the large deviation functional (31) is 
quadratic [32,35]. In general the complicated expression of the cumulants (23) expresses the non-trivial coupling 
between the flux j and the optimal density profile po(x). 



4.2. The cumulants 



The parametric expression (34, 35) for the large deviation function Fl{]) can be transformed into another 
parametric form for Pl(X) defined in (4,5) 



p, L (\Pa,Pb) = ~ 



f D( 

J x/T+ 



D(p) dp 



y/l + 2Ko{p) 



(37) 



with K = K(\, p a , pt) is the solution of 



-K 



DM 
(p) 



v/1 + 2Ka(p) 



(38) 



By eliminating K, one can obtain the expansion of fj,L in powers of A and by taking successive derivatives (5,6) 
with respect to A one gets for the cumulants of the current: 



(Qt) i r 

— = L h > 



(Qi) - {Qt? 



(Q'i)c 



Ik 

t LIi' 
1 3(7 3 7i - if) (Qi) c _ 1 3(5/ 4 J 1 2 - Uhhh + 9/|) 
" L 



(39) 



t L If ' t 
where the integrals /„ are given by 



Pa 

I n = j D(p) a(p)"- 1 



dp . 



In the case of the SSEP, one has D(p) = 1 and a(p) — 2p(l — p). One can simplify (37,38) and get 



Mi (A) = j [\og(VTT^ + V^)] with uj = (e A - I) p a + (e - l)p b - (e x - l)(e - l)p a p b - 



(40) 



From this one can recover the cumulants (23) and determine all the higher cumulants. Expressions equivalent to 
(40) were derived in the theory of shot noise of mesoscopic conductors [36,37]. 

Remark. The density p is the physical relevant parameter. However it can be useful [38] to consider instead the 
conjugate field f3 = log 2 (see (21)). Formula (31) then simplifies as it depends only on one macroscopic input cr(/3) 



F L {j,Pa,Pb) = min y / 





dx 



2ff (/?(*)) 

where the minimum is taken over the chemical potential profiles such that /3(0) = logz a and ,0(1) = log 25. 



(41) 



4.3. Heat flux. 



All the above discussion can be generalized to the case of a heat flux in diffusive systems: one has to replace 
everywhere the density profile p(x) by the temperature profile T(x). There is even one simplification as in the 
thermal case D and o are related as in (17) so that for D(T) defined as in (30) (i.e. (Qt) /t = (T a - T b )D(T)/L 
for small T a — T b ), one gets 
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Lj 



-t a 

-I 



dTD(T) 



[1 + 4KT 2 D{T)] 1 / 2 



_ . f dT 
- 3 J 2T 2 



1 + 2KT 2 D(T) 
[1 + AKT 2 D{T)] 1 / 2 



- 1 



(42) 



5. The macroscopic fluctuation theory 



Building on the hydrodynamic large deviation theory [39,29,30], Bertini et al developed [40,41,42] a general 
framework to determine the steady state large deviation function of non equilibrium systems. This framework has 
been extended [12,43] to the current large deviations. Let us sketch briefly their approach. For diffusive systems 
(such as SSEP), the total flux Qi(t) flowing through position i between time and time t and the density pi(t) 
near position i are, for a large system of size L and for times of order L 2 , scaling functions of the form 



Qt(t) 



and 



Pi{t) 



p ^vT 2 



It is convenient to introduce the instantaneous current defined in terms of the rescaled time r 



q(x,r) 



dQ(x,r) 
dr 



(43) 



In fact Lq(x, r)dr is simply the total flux of particles through position [xL] during the microscopic time interval 
[L 2 r, L 2 (r+dr)], with 1/L <g; dr <g; 1 so that there is a large number of particles which contribute to the integrated 
current but the density does not vary over this small time interval. Remark that the current q is defined after a 
diffusive rescaling, i.e. the space is scaled by 1/L and the time by 1/L 2 . Thus unlike the microscopic current, q 
remains of order 1. The conservation of the number of particles implies that 



dp(x,r) _ d 2 Q(x,r) 



dr 



drdx 



dq(x,r) 
dx 



(44) 



The macroscopic fluctuation theory [12,43] gives for the probability of observing a certain density profile p (x, r) 
and a current q (x, r) over the rescaled time interval < r' < r 



Pro^{p(a;,r'),(7(a;,r')}^ ~ exp —L J dr' J 



I } \g(x,r') + D( P (xy))^lV 
' ■' - ' dx- 



2a(p(x,r')) 



(45) 



Of course p and q have to satisfy the relation (44). (Note that if t is the microscopic time, then r = t/L plays 
the role of a macroscopic time). A similar expression was obtained in [46,47] by considering stochastic models in 
the context of shot noise in mesoscopic quantum conductors. The functional (45) was used to calculate the large 
deviation functional of the density for several systems [41,44] and in the case of SSEP the results agree with an 
exact microscopic derivation [45]. 



The large deviation function (2) for observing the total current j, i.e. the following event 

T 1 

Lj — — J dr' J dx q(x, r') 
o o 

as predicted in [12,43] by the macroscopic fluctuation theory (45) becomes 

[q(x,r') + D(p(x,r'))^lf 



(46) 



J-(j) = — lim — min 

L T^OO T p{x,T r ) 



dr' J d.r- 

o 



2a{ P (x,T>)) 



(47) 



where the minimun is over all the density profiles {p(x,t'),0 < t' < t} and the current {q(x,T'),0 < r' < r} 
which satisfy the conservation law (44) and the global constraint (46). 

If the optimal density and current profiles are time independent (up to boundary effects for r' close to or r 
which do not contribute in the r — > oo limit), one recovers the predictions of the additivity principle (31) and 



11 



= When the optimal profile is time dependent the additivity principle predictions (31) give only an 

upper bound : !F(j) < F L (j). In [12,43], Bertini et al provided an example for which the functions F L and T are 
different. In their example, Fl was not a convex function of j and T was its convex envelope (see (3)). They also 
proved (see [43] Section 6.1) that T reduces to Fl under the following global condition on D and a 

For all p, D(p)a"(p) < D'(p)a'(p) (48) 

This holds for the SSEP and the Zero Range process [30]. Condition (48) is however only sufficient. 

In general a dynamical phase transition [43,48] may occur where the system switches from a time independent to 
a time dependent optimal profile. To calculate the large deviation function T one needs to determine the optimal 
time dependent profile p(x,r), which is not an easy task as the optimization problem is non-linear. A complete 
characterization of the regime for which the additivity principle holds (T(j) = Fl(3)) remains a challenging 
problem. 



6. Phase transitions 



In this section, we try to determine the phase boundary where the optimal profile becomes time dependent. To 
do so we consider a more general situation with a small driving force (like an electric field) in the bulk. Adding 
such a driving force of amplitude v/L to an open system of length L with reservoirs p a ,pb modifies the mean 
current (30) as follows [29,48] 

^ L ^{Pa-Pb)^ L + ^a(Pa) for p a - p b small (49) 

In (49), the conductivity a which was defined as the variance of the current in (29) can also be understood as the 
linear response to the small field v/L. The effect of the field can be easily taken into account in the framework 
of the macroscopic fluctuation theory [43,48] by arguing that locally the current has Gaussian fluctuations with 
mean value given by (49). The functional (47) becomes 



[q(x, t') + D(p(x, T '))22^2 - va(p(x, r'))] 



2 



Hj) = -r lim - min / dr' / dx± -—- -r ] — . (50) 

J ' L t >oo t p{x,t') j J 2<r(p(x,T')) v ; 



40»,t') 

with the constraint (46) on the total current. 



6.1. Stability of the functional 

The time independent optimal profile po(x) is now a solution of (see (32)) 
(D{p {x))p' {x)) 2 = (jL - va(p (x))Y + 2Ka(p (x)) , 



(51) 



where the constant K has to be adjusted so that po satisfies the boundary conditions (po(0) = p a and po(l) = Pb)- 
A situation for which po(x) is certainly not optimal in (50) is when a small time dependent perturbation is 
sufficient to lower (50). To investigate the stability of po(x) against such perturbations, one can write 

p(x, t) = po(x) + Sp(x, t) 

q(x,r') =j + Sj{x,r') 

where Sp and Sj have zero time averages and are related by (44). Inserting these expressions into (50), one gets 
at the second order 





where we introduced the functions 



A'{po)p'o 



3pSp' + - 



B"( P0 )+A"( P0 )(p' ) 2 



(8 P y 
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A(u) 



2a(u) ' 



and B(u) 



2^0 + ^ a(U) 



The coefficients of the quadratic form are x dependent but time independent. Thus to analyze the stability, one 
can consider perturbations of the form 



— exp(iiL>r')ipi(x) — — cxp(—iLL>T')ip' 2 (x) 

LJ LJ 



5p{x,r') 

5j(x,T r ) = exp(iLJT')tpi(x) + exp(—iujT')ip2{x) 
where tp2(x) = <P*(x). The quadratic form can be rewritten as 
1 

f , (<fii<fi2 . cr'(po) i r / n . 2 „ „ 

J { v{Po) cr 2 (po) OJ L OJ 2 



(52) 



(53) 



^'(po)po 



r ' "1 "Mi 
|_¥>i <P2 + <Pi <Pi\ + 



B"(p )+A"(p )(p' y 



For the time independent profile po(x) to be stable against small time dependent perturbations, the quadratic 
form (53) has to be positive for all lj. In general, all the coefficients in (53) are spatially dependent through po{x) 
and it is difficult to provide from (53) a more explicit characterization of this local stability. 



6.2. Periodic systems 



For a system of N particles on a ring of length L with density p = N/L, the flat profile po(x) = p remains a 
solution of the Euler-Lagrange equation associated to (50). In this case, the coefficients of the quadratic form have 
no x dependence and the different spatial modes decouple. One can choose ifii = exp(i kx) and if2 = exp(— i kx) 
and the positivity of the quadratic form implies that for any k (multiple of 27r) and lj, one has 



1 



1 + Lj 



ko'jp) 
Loa(p) 



+ — 



k 2 Djpf 



a"{P) , 2<r"(p) 

MP) 2 



+ V 



> 0. 



(54) 



The first mode to become unstable is the fundamental mode k = 2n, thus the flat profile is stable when [48] 



8tt 2 D 2 (p) 
°{P) 



> v"(p) 



(Lj) 2 



Y° 2 (P) 



(55) 



Let j c be the critical current for which (55) becomes an equality. If o~"(p) > (resp cr"(p) < 0), the flat profile 
becomes unstable for currents \j\ > \j c \ (resp \j\ < \j c \)- Remark that the instability regime is always symmetric 
with respect to as predicted by the fluctuation theorem (19) which becomes in presence of a driving force 



Hj) - H-j) = -2vj -2j J ° 



Dip) 
(P) 



dp . 



Beyond the threshold (55) a bifurcation occurs and a traveling wave of the form p(x — vt) is more favorable than 
the flat profile p. From (54), we get that close to the phase transition, the optimal velocity is given by 



lj — —2nLj 



.. °'(P) 



°(P) 



v = Lj, 



°'(P) 
o(P) 



(56) 



If we make the assumption that no first order transition occurs before the second order transition predicted at 
j c , we can compute the expansion of T close to j c . We consider small current perturbations 

Lj = Lj c + e, with Lj c = ^jvMP) 2 + 8 ^iy ^ (57) 

Let us limit the discussion to the case cr"(p) < 0. Then for e > the flat profile remains optimal and one expects 
that the large deviation function is quadratic 
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Ve >o, m = w < + lm" m * < 58 » 

On the other hand for a"(p) < 0, the flat profile is unstable for e < and the expansion of !F(j) at the second 
order in e can be obtained by approximating the travelling wave as follows 

po(x,t) w p + \/— eai sin (2-k(x — vt)^j + ea 2 cos ^4tt(x — vt)^j (59) 

where the velocity v = Lj c a a (S is given by (56). In fact the expansion should also include corrections in e to the 
velocity as well as other Fourier modes, but a computation shows that they do not contribute to the second order 
expansion of T. Inserting the test function (59) in (50) implies that the second order of T(j) in e is given by the 
quartic form 



Ljco"{p) , / 3ir 2 D'(p)D(p) , 3ir 2 D(p) 2 a'(p) , tt 2 D{pfo^ (p) 



2 



i*(pf + \ oift) + 2a(p) 2 + 2a(p)a"(p) 1 " J 1 " ] 

n 2 D(p) 2 a^(p) n 2 D(p) 2 a'(p) 2 n 2 D(pf a"(p) K 2 D(p)D'(p) a'(p) 



16 a(p)a"(p) ' 4a(p) 3 4a(p) 2 2a(p) 2 

~ 2 D(p)D"(p) 7v 2 D'(p) 2 3v 2 *"(pf 

4a(p) + 4r(p) + 64a(p) y ' <r(p) ^ ' 2a(p) 

where a^ z \a^ denote the third and fourth derivatives. The optimal amplitudes a\,a 2 of the traveling wave (59) 
are the minimizers of the quartic form. Note that (60) is not always stable and for some specific choices of the 
functions D and a, the mimimum of (60) can be -co. For example the sign of the coefficient a\ depends on <j^(p) 
which can a priori take any arbitrary value. The condition for the quartic form (60) to be stable can be written as 



>0, (61) 
where 

A(p) = 9u 2 a(p) 2 a"(p) 4 - 96n 2 D'(pf a (p) 2 a" (pf + 48tt 2 D(p)D" (p) a (p) 2 a" \pf 
+ 127T 2 ■D{p) 2 o'{pfo"{pf - 24iv 2 D(p) 2 a(p)a'(p)a"(p)a^(p) 

+48n 2 D(p) 2 a(p)a"(pf + 48n 2 D(p)D' (p) a (p) a' (p)a" (p) 2 + 48n 2 D(p)D' (p)a(p) 2 a" (p)a ( - 3 \p) 
-^ 2 D{p) 2 a(p) 2 o^\p) 2 - 127V 2 D{pfo{p) 2 o"{p)cj {4 \p) . (62) 
If (61) is not satisfied then one expects that a first order transition occured before j c . 

We suppose as before that cr"(p) < and that (61) is satisfied. Then the minimum of (60) is achieved for 



ai= 6L jM p)_ : "m 3 and a^(m^M-^k)^ (63) 



V2D(p) 4cr(p) 12(7" (p) 
where A(p) is defined in (62). Thus for e < 0, one finds at the second order 

m = (Ljc + £ ;r ( ' ))2 - 3a y^ J ?c 2 +oU) . (64) 

2a(p) <t(p)Hp) V J 

Comparing (58) to (64), we see that if A(p) > the time dependent profile gives lower than the flat profile. 

As an example, we consider the weakly asymmetric simple exclusion process (WASEP) which follows the same 
exclusion rule as the SSEP with jump rates biased by exp(^/L) to the right and exp(— v/Li) to the left. In this 
case D(p) — 1 and er(p) = 2p(l — p). Since cr"(p) = —4, the threshold of stability of the flat profile is given by 

Lj c = ^Jv 2 o{p) 2 - 2tt 2 ct(p) 
For e < 0, a perturbation of the form (59) leads to the quartic form (60) 
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(65) 



The quartic form (65) is stable so that the amplitudes of the traveling wave (59) are given by 



2Lj c a(p) a'(p) 2 

" 1 ^/ > Mp?-6^(p) + ^ and a2 = ~MP) ai (66) 



Finally for e < 0, one finds at the second order 

(Lj e + £ -^(p)) 2 _ ^ 2 q(p) - 2^ 2 a / 3 N 

KJ > 2a(p) 3is 2 a(p) 2 -6it 2 tr(p)+n 2 + V ) 



2a(p) 3v 2 a(p) 2 — 6n 2 a(p) + n 2 

It is interesting to note that for the WASEP on a ring at density p = 1/2, the optimal velocity is 0, thus the 
optimal profile remains time independent. This means that the phase transition could have been detected already 
at the level of the functional Fl ■ 

Finally, let us mention that in the large drift limit v — > oo, the asymptotic cost for (50) as well as the asymptotic 
shape of the optimal traveling waves can be computed [48]. In particular for the WASEP, the current large 
deviation function (50) converges in the large drift limit v — > oo to the current large deviation function of the 
totally asymmetric simple exclusion process [49,50]. We refer the reader to [48] for further details and to [51] for 
a study of the large drift limit in the case of open systems. 



7. Conclusion 



In this paper, we have shown how to generalize the detailed balance relation to take into account the effect of 
reservoirs. From this generalized detailed balance relation the fluctuation theorem [16,17] can be recovered which 
characterizes the odd part of the large deviation function of the current. By a simple additivity principle [35], one 
can predict for diffusive systems the whole large deviation function as well as all the cumulants of the current. 
These predictions agree with previous exact computations for some stochastic models like the symmetric simple 
exclusion process [32]. For some models, however, the additivity principle provides only an upper bound of the 
large deviation function of the current [12]. This fact as well as the occurrence of phase transitions has been 
discussed in the framework of the macroscopic fluctuation theory [43,48]. 

A challenging issue would be to characterize precisely the range of validity of the additivity principle in the 
case of diffusive stochastic models. Here we were only able to address the local stability of the time independent 
solution of the additivity principle. How to calculate the large deviation function of the heat or particle current in 
a more general framework (several species of particles, additional conserved or non conserved quantities) is also 
an interesting open question. 
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